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Abstract 

We investigate the interactions of collective excitations with vor- 
tices in superfluid systems, including ^He and superconductors. The 
dynamical equations are obtained by the aid of the many-body wave- 
function and the density-density correlation function. The scattering 
cross-section of collective excitations with a vortex is calculated in the 
Born approximation (valid at long wavelengths), and is expressed in 
terms of the collective excitation spectrum a;(g) by the simple and 
general formula a = f ■p/^p' cot^ |, where q is the wave number of 
the excitation, and 6 is the scattering angle. At short wave lengths, 
the classical equations of motion are derived. 
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1 Introduction 



The motion of vortices is considered to be the main source of dissipation in 
superfluid systems. At finite temperatures, the existence of collective exci- 
tations and their interactions with the vortex lines create the phenomenon 
known as mutual friction. |], Q The phenomenological parameters of 
mutual friction can be derived from the 2D scattering cross-section of the 
collective excitations P). 



Different methods have been used in the literature for the calculation of 
the scattering cross-sections. One method is based on the classical dynam- 
ics of phonons or rotons ^, |^ treated as particles. In the presence 
of a vortex, the Hamiltonian for such a particle is given, for instance, by 
H = eo(p) + V{r) ■ p where eo(p) is the energy spectrum of the collective 
excitations, and the second term represents the Doppler shift of energy due 
to the velocity field V{r) associated with the vortex. Effects of the vortex 
core and density modulations have also been considered. Another method 
is to treat such a particle quantum mechanically by introducing a quantum 
wave function for it, using the operator version of the above Hamiltonian 
to describe the dynamics. This method has been applied to the roton case, 
with the cross-section calculated in the Born approximation There is 
then the more conventional method based on the hydrodynamic wave equa- 
tions. Born approximation ^ and phase shift analysis have been used 
for the calculation of phonon scattering cross-sections. 



We observed that the phenomenon can be considered in more general 
terms. Besides ''He, there are other systems, such as superconductors, 
quantum Hall systems , quantum spin systems, that support vortices and 
collective excitations with different microscopic structure and ground state 
properties. The purpose of this paper is to establish the relation between the 
different methods, and to present a calculation of the scattering cross-section 
that does not refer to the details of the microscopic structure of the systems. 
For long wavelength scatterings, the cross-section may be calculated using 
the Born approximation, and is found to be given by the simple and general 
expression: 
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where q is the momentum, uj is the excitation spectrum, and 9 is the scat- 
tering angle of the excitations. We also show that a classical interpretation 
of the scattering is possible in the opposite limit, i.e. at short wavelengths. 



In Section 2, we will give the basic ingredients of our method. In Section 
3, we will calculate the cross-section in the Born approximation, and present 
the classical approach. In Section 4, we will apply the results we have derived 
on some systems. In Section 5, we present our conclusions. 



2 Dynamical Equations 

In this section we will lay out the general formulation of the problem, and 
obtain the dynamical equations for the elementary excitations in the presence 
of a vortex. The core of the idea is to obtain an Hamiltonian for the elemen- 
tary excitations, and use the Hamilton's equations to describe the motion of 
the excitations. 



2.1 General Formulation 

The elementary excitations are represented by small variations in the lo- 
cal density and phase, and are described by the following many-body wave 
function [Tl 



/ N 

^' = exp ^ 
\i=i 



a{xi) + iS{xi) \E'o, 



(2) 



where \&o is the ground state wave function, S and a are arbitrary func- 
tions that satisfy a{x) -C 1, and do not vary rapidly over the characteristic 
length scales of the system, such as the interparticle spacing in ^He and the 
coherence length in superconductors. 



We follow the lines in Jill to obtain the energy in this state asQ: 



E-E, = -lp{x){ WS{x) 



Va(f) 



dx. 



(3) 



^we use % = m ^ 1 unless otherwise stated 
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where the integration is over the spatial dimensions of the system, and p{x) 
is the density: 



p(f) = (^lp(f)l^) 

= (^o| p(x) exp (^2 J p{x')a{x)dx^ l^^o) , (4) 

with p(x) = J2i ^{^ ~ ^i)- 

We want to express eq. @ in terms of the canonically conjugate variables 
6p = p — Pq and S JTSf, and regard it as the Hamiltonian of the collective 
excitations. Eq. (|^) is exact, but we have to make an approximation in 
order to make progress. Thus we will linearize the density in a using the 
assumption that the fluctuations are small, yielding 

p{x) = po + 2 1 {^o\p{x)p{x')\^o)a{x')dx (5) 
= po + 2po / g{x — x')a{x')dx' , (6) 



where po is the density of the system in the ground state, and the function 
g{x — x') is the density-density correlation |T^]. We can invert this relation 
and solve for a{x) 

^i^) — o — / ^p{^')9^^{^ ~ x')dx' , (7) 
with the definition of the inverse function by 

g'^{x - x')g{x')dx' = 5(f), (8) 



Relation (|^) enables us to write the energy function eq. (^) in terms of 
the variables Sp{x), S{x) and the correlation function g{x). The properties 
of the ground state and the interactions between particles are included in 
g{x). The normal modes of these equations will be the excitation spectrum 
of the fluctuations. 
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2.2 Collective Excitations in the Absence of a Vortex 



After substituting eq. (0), the Hamiltonian of the system without any vortex 
in it can be written as 



H 



+ 2^0 



VS{x] 



dx 



1 

Wo 



5p{x')g ^{x — x')dx' 



dx. 



(9) 



The first term can be considered as the kinetic energy and the second one 
as the potential energy of fluctuations. The Hamilton's equations are just 
functional derivatives of the Hamiltonian : 



d Sp{x, t) dH 

Tt ~ dS{x, t) 

= -poV^Six,t) (10) 

dS{x,t) dH 



dt dSp{x,t) 

= — ^V' / f g-\x-x')g^\x' -x")6p{x")dx' dx",{ll) 

where we have used the fact that for an isotropic system g{x) depends only 
on |x|, and omitted some boundary terms. If we recall that VS is the velocity 
field, we can easily recognize the first equation as the continuity equation, 
and the second one as the generalized Euler equation of hydrodynamics in a 
form that also takes into account of density correlations in the system. 

We can try a plane wave solution to these equations 

S{x,t) = Ae*(^-^^-"*) 

Spix,t) = 5e'(^-"^-"*), (12) 
which immediately requires the conditions 

-iBuj = Apok"^ (13) 
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iAu 



B- 



4po 



n 2 



B 



(14) 



where S{k) is the structure factor defined as the Fourier transform of g{x). 
A non trivial solution to these conditions is 



(15) 



with 



B 



2ipoS{k)A. (16) 

The first relation is the famous Bijl-Feynman formula for the collective ex- 
citation spectrum []T2[, and the second gives the relation between the phase 
and amplitude fiuctuations. 



2.3 Effective Hamiltonian in the Presence of a Vortex 

Thus we have shown that eq. (^), when used as the Hamiltonian of the 
collective excitations, gives us the correct spectrum. Now we want to put 
a vortex into the system and obtain an effective Hamiltonian for the long 
wavelength excitations. For simplicity we assume a single, straight, fixed 
vortex line at the origin. The many body wave function in eq.(l) can now 
be written as: 

/ N \ N 

VI/ = exp a(f,) + tSix,) n / (l^. l) e^'^^o, (17) 

\i=l / j=l 

where / is a function that becomes 1 outside the vortex core but goes to 
zero at the origin. 6j is the polar angle at the position of particle j. A full 
treatment of the problem should take into account the core structure through 
the function /. However we will make the further simplification of taking / 
equal to 1 everywhere, thus disregard the effects of the vortex core. Then the 
effect of the vortex is to simply cause a shift in the canonical phase variable 
by^ 

VS^VS + -, (18) 
r 
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where 6 is the polar angle and r is the distance to the vortex line. Under 
these circumstances, the Hamiltonian up to second order terms in 6p{x) and 
VS{x) is 



H ^ - 




yzpo J 



where we have ignored some boundary terms and terms second order in 
the vortex field. As discussed above, the first and last terms describe the 
collective excitations. The scattering of excitations is caused by the second 
term which, in classical terms, looks like a Doppler shift in the energy. The 
correlation function g{x) is unchanged by the inclusion of the vortex as the 
function / representing the vortex core was taken to be 1 everywhere. For 
the same reason the vortex couples to the kinetic term but not the potential 
energy term. This is why the vortex interacts with the excitations in quite 
a universal way, as a result of which a simple and general expression for the 
scattering cross-section will be obtained in the next chapter. 



3 Calculation of the Scattering Amplitude 



Now we can analyze Hamilton's equations, and calculate the scattering am- 
plitudes. The Hamilton's equations in this case differ from the previous ones 
only by the appearance of the vortex terms 



The time dependence of the solutions can still be factored by e then 



dSp{x, t) 
di 



-poV^S{x, t)--- V5p{x, t) (20) 
9 - 

-■VS(x,t) 
r 

- -J-^^ I I ^"'(^ - ^)^"'(^' - x")5p{x")dx' df".(21) 



dS{x, t) 



dt 
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the equations read 
— iu}5p{x) = 



e 



(22) 



+ — V / / g-\x - x')g-\x' - x")6p{x")dx' dx".{23) 

The vortex terms have position dependent coefficients and prevent an 
exact solution to be obtained. We will use a Born approximation, in which 
we will assume that we can replace the S{x) and 6p{x) in the vortex terms 
by the free solutions obtained previously: 



So{x) = Ae"^-'^ 
Spo{x) = Se*«-^ 



(24) 



with the same dispersion relation. We will also restrict ourselves to two 
dimensions in the plane perpendicular to the vortex line, since there is no 
scattering parallel to the line. The dynamical equations now contain inho- 
mogeneous terms 



iuSpix) + poV^S{x) ^-iB-- qe'^-^, 

r 



(25) 



iu!S{x) + 



4po 







g-'(x-x')g-'(x'-x")6p(x")dx'dx" = ^--ge^^'^ (26) 

r 



By taking Fourier transform of both sides we turn the coupled equations into 
a linear matrix equation 



-k'^Po 



452(A;)po 



lUJ 



m 

5p(k) 



r 



A 



B 



(27) 



where a tilde denotes the Fourier transformed quantities. It is simple alge- 
bra to perform the necessary Fourier integral on the right hand side of the 
equation. The result is 



2tt 



z ■ {qx k) 



q — k 



2 ' 



(28) 
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where z is the unit vector in the direction of the circulation. We invert 
eq. (pTl) , replace iS(g) by g^/[2ci;(g)], and then obtain the Fourier transform 
of the scattering amplitudes 



Sp{k) 



A2 



TT 



z ■ {qx k) 



-LJ^{k)+u^{q) 



q — k 



luo 



Pok"^ 



pok'^ 



(29) 



The transformation back to normal space variables is given by 



d^k 



- S{x) - 


= J e^^-^ 


" S{k) ' 






6p{k) 



(30) 



The integrals cannot be evaluated exactly. Since our purpose is to calculate 
the cross-section we only need the amplitudes in the radiation zone, /cr 3> 1. 
A separation of the k integral over components, k\\ and k±, parallel and 
perpendicular components of to a; respectively, is helpful. Then specifically 
for S{x) we get: 



S{x,t) 



Aie 



27r 



dkw dk I 



^iknr 



(27r)2y uj{k)-uj{q)uj{k)+uj{q) 
uj'^{k)q^ + uj'^{q)k'^ z ■ {q x k) 



uj{q)k'^ 



k 



(31) 



We will first keep k± constant, and evaluate the k\\ integration. We can 
take care of the boundary conditions at infinity by supplying a small imag- 
inary part wherever necessary. The boundary conditions require that the 
outgoing scattering amplitudes should behave as e^'^^/i/r, so the relevant 
poles are at the values of k\\Q = K"^ — k'j_ + ii], where K satisfies 

u;{q) = u;{K), (32) 
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and ?7 is a positive infinitesimal constant. 

We have to note that k is not necessarily equal to q, as might be expected 
from energy conservation, since the spectrum may exhibit a behavior such 
that several k values can correspond to the same energy. A very well known 
example is the roton spectrum, which we will investigate in the next section. 
For the time being we will work with the case K = q. 

Expanding around the pole to linear order in {k^\ — kp) gives 

dk 

uj{k) ^ u{q) + ^'(g)^ fc,|o {h - ho)- (33) 

Then we place it into the integral and obtain the simple result 

= -Ae-^"*Res(q), (34) 



with 



Res = - / dk^e'^W'^' ^— 



g(fc|[osin6' - k^ cos 9) 

^ kj^ + kl + q^ -2qk\iocose -2qk^sme' ^ ^ 



and keeping in mind that fcyo = \J q^ — k\. Now we have to evaluate the 
integral over k±. In the radiation zone where we take r — > oo, the integral 
will be dominated by the stationary point of the exponential, which is lo- 
cated at fc^ = 0. Assuming that the integrand other than the exponential 
behaves smoothly, we can set fc^ = everywhere except in the exponential. 
A stationary phase approximation gives 

^ ,-7r nre^^^ ^fq sin^ , „. 

Res = -e-*4./ -^{^ -. (36) 

V 2 v^cu'(g)l-cos^ ^ ' 



Combining with eq. (0) we get the desired result 



^ ' ' V 2uj\q) 1 -cos^ ^ ' 
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The cross-section is given by the square of the amphtude of the outgoing 
wave 

^(9,^^) = f-fT^cot^|, (38) 
2 [uj'{q)] 2 

which is the central result of this article. Note that the angular dependence 
is universal, which is symmetrical with respect to ^ — —6, and diverges at 
small angles. The spectrum function enters in a simple form, and only affects 
the q dependent factor. 

One can evaluate the integrals in the general case where there may be more 
than one possible value of K. A similar calculation leads to the scattering 



amplitude given by the sum of residues satisfying eq. (|3^) 



e'^'' nrq{K^ + q^) 1 sin 9 



K V 2 ^ uj'{K)K^ + q^-2Kqcos9' 

(39) 

to pay more attention to the degenerate case when uj'{K) = 0. In that 
case the idea is still the same, only we need to expand to quadratic order 
in eq. (0), which forces us to evaluate the derivative of the integrand in 
eq. dH). 



Born approximation is valid if the scattered amplitude is small when com- 
pared with the incoming amplitude. We can expect this to be violated most 
where the potential is strongest, which is the origin in our case. A simple 
calculation after setting x = in eq. (PT|) leads to the requirement that 

^ « 1. (40) 



uj'iq) 

Thus if the excitation spectrum vanishes with a power less than 2 as g 
then Born approximation is valid in this limit. 

3.1 Classical Hamiltonian 

The other approach which is mostly used for rotons is the classical equations 
of motion for the collective excitations. It is possible to show that one can 
obtain the classical Hamiltonian, thus the classical equations, from eqs.(pO|) 
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and (0). These equations are like wave equations for the amphtudes S and 
6p. The classical equations follow from the geometrical optics approximation 
or the WKB. We replace the amplitudes S and 6p in the eqs.(^D|) and (PH) 
by 



S{x, t) 
Sp{x, t) 



(41) 



where A and B in general have spatial variation, but in the spirit of the 



WKB approximation they can be taken as constant. Then from eqs. (|20|) and 
ni) we obtain 



—I 



+ - ■ V0 

at r 



+ - ■ V0 

ot r 



B 



A 



Po 



V0 



V0 



A 



4po52(V0) 



-B. 



(42) 



(43) 



We want to remind that we are stil using units h = m = 1. The momen- 
tum and energy of classical particles corresponding to wave packets of the 
amplitudes are given by 



E 



Thus we obtain the relation 



90 
"dt 



and p = V0. 



(44) 



E = uo{p) + - ■ p, 



(45) 



where we identified the energy spectrum of the excitations, u!{p) as in the 
Bijl-Feynman formula. This is obviously the classical Hamiltonian for the 
collective excitations with the Doppler shift in the energy due to the vortex. 
The equations of motion that follow from this Hamiltonian have been used 
extensively in the literature for rotons [0, 0, ^, and recently for scattering of 
phonons from a classical vortex 0. 

The WKB approximation is valid, if the spatial variation of t) satisfies 

V^<p < |V0|2, (46) 
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For this approximation to be valid the wave packets have to satisfy 

V-p«|p1% (47) 



or as an estimate on the left hand side we can use the uncertainties of mo- 
mentum and position 

(48) 

We can take Ap ~ 1/Ax and Ax ~ C,- The classical approximation is valid 
if 

(49) 

Thus WKB is apropriate for large momentum. This is the opposite limit of 
validity of the Born approximation given in eq. (|40|) . 



4 Applications 

The results of the previous section are quite general in the sense that as long 
as we know what the collective excitation spectrum of a system is, we can 
find the scattering cross-section of these excitations from a vortex. In this 
section we will show their use in superfluid helium and superconductor thin 
films. 



4.1 Superfluid He 

1. Phonon Scattering : In the phonon part of the spectrum, the long wave- 
length low energy part, the relation between k and lj is one to one, and is 
given by 

uj{q) = sq, (50) 

where s is the speed of sound. The only pole satisfying uj{q) = uj{K) is at 
K = q. Thus from eq. (|38|) we obtain: 




^(9'^)=o? — cot^o- (51) 



This result was obtained by Fetter using phase shift analysis of the 
scattered wave in the hydrodynamic equations. The other results in the 
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literature vary. Pitaevskii[^] and Sonin0 have different results from the 
above although they both use Born approximation and similar equations. 
However we have found that Pitaevskii's calculation suffers from an algebraic 
error. In his paper eq.(25) does not follow from the substitution of eq.(24) 
into eq.(19). Instead when correctly done the result is the same as what we 
have obtained above. The calculation of Sonin is a more complicated one 
that it takes into account of the curvature of the vortex line. His calculation 
must be also incorrect, however, as it gives the result of Pittaevskii in the 
limit of a straight vortex line. 

2. Roton Scattering : We can employ the formula (^) to obtain the cross- 
section of rotons. There are two values that satisfy eq. (P^. The spectrum 
near the roton minimum can be written in the Landau form: 

^,) = A+<i^, (52) 

where A denotes the energy and go denotes the wave vector at the roton 
minimum. The roton has two poles that satisfy w{q) = u){K), which are 

Ki = q and K2 = 2qo — q. (53) 

Now using eq.(^) : 

Six,t) = Ae-'Jl'^-^n 
\ 2 y/r \q- qo\ 

.J sing sing ] 

\l-cosg 2go-gi-^giH|^cosg/- ^ 

We immediately notice that the scattered wave has two parts, one spherical 
wave with momentum q and one echo with 2go — q- The cross-section will not 
only involve the individual contributions of each part but also the interfer- 
ence. For |g — gol ^ ~ the interference can be ignored and the cross-section 
will be given by: 

This is also the same result obtained long ago by Hall and Vinen0. 
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4.2 Superconductor Thin Films 



Another system where we can use formula (^) is a superconductor thin film. 
Scattering cross-section and mutual friction in superconductors were investi- 
gated on the basis of quasi particle-vortex interactions. This is the dominant 
contribution for ordinary superconductors. However, at low temperatures 
and small core size, the quasi particles are frozen, and interactions between 
vortices and collective excitations become more important. 



In a very broad sense, superconductivity is superfiuidity of charged parti- 
cles. The important distinction is that the vortex singularities are like thin 
fiux tubes for bulk materials because of the diamagnetic currents set up in 
the sample. However, in thin films this screening is only effective beyond a 
length Lg given by : 

i.^f (56) 

where Aj;, is the London penetration depth and d is the thickness of the film. 
This length scale can be quite large compared to the bulk value A^,. Thus 
the vortex in a thin film is very much like a vortex in neutral ^He. There is 



one little difference that is related to pairing in the ground state [0]. This 
forces us to introduce a modified shift of ^ in equation ([181) . The final effect 



of this on the cross-section is a factor of \: 

The collective excitations of a superconductor are the density fluctuations 
of cooper pairs, and have the same dispersion as the ordinary plasmon at zero 
temperature, which has the following form in two dimensions: 



u;{q) = (58) 
\ ) 

where n is the 2D density of electrons, e the electron charge and me the 
electron mass. The resulting cross-section is then: 

a{q,e)=^(^\ L,q\ot^^-, (59) 
where c is the speed of light. 
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5 Conclusions 



In summary, we have investigated the interaction between collective exci- 
tations and vortices in superfiuid systems. Starting from a Feynman type 
many-body wave function, general dynamic equations for the collective ex- 
citations in the presence of a vortex have been derived, in which the spe- 
cific properties of a given system enters only through the density correlation 
function or the static structure factor. The scattering cross-section has been 
derived in the Born approximation for long wavelengths, whose angular de- 
pendence has the universal form cot^ |, with a prefactor only involving the 
energy spectrum of the collective excitations. For short wave lengths, the 
geometrical optics or the WKB approximation to the dynamical equations 
has been shown to yield the classical equations of motion for wave packets. 

We have applied the Born approximation cross-section formula to phonons 
and rotons in ^He, and showed that the results are consistent with previous 
calculations. As a new result, we have also obtained the cross-section of 
plasmons from vortices in superconductor thin films. These cross-sections, 
when combined with the formulae in can be used to find the frictional 
forces on a vortex. For the case of superconductor thin films, the longitudinal 
dissipative force on a moving vortex is shown to be proportional to a high 
power of temperature, T^°, which supports the idea that vortices move freely 
at low temperatures. 

The transverse force due to scattering vanishes as long as the Born ap- 
proximation is valid, which is the case for phonon and plasma scatterings. 
However, the dynamics of rotons in superfiuid helium falls rather into the 
validity range of WKB approximation, and can yield a transverse force on 
vortices ^, |^. This is very similar to the scattering of electrons by a mag- 
netic flux line. The classical calculation valid at short wavelengths yields 
a transverse force, whereas the Born approximation valid in the long wave- 
length limit does not. 
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